Remote Sensing (RS)
Hosein Nesari; Reza Shah-Hosseini; Amirreza Goodarzi; Soheil Sobhan Ardakani; Saeed Farzaneh
Abstract
Extended Abstract
Introduction
Atmospheric aerosols are a colloid of solid particles or liquid droplets suspended in the atmosphere. Their diameter is between 10-2 to 10-3 micrometers. They directly and indirectly affect the global climate by absorbing and scattering solar radiation, and they also ...
Read More
Extended Abstract
Introduction
Atmospheric aerosols are a colloid of solid particles or liquid droplets suspended in the atmosphere. Their diameter is between 10-2 to 10-3 micrometers. They directly and indirectly affect the global climate by absorbing and scattering solar radiation, and they also have a serious impact on human health by emitting harmful substances. In addition, high concentrations of aerosols on a local scale due to natural or human activities have adverse effects on human health, including cancers, pulmonary inflammation, and cardiopulmonary mortality. Monitoring the temporal and spatial variability of high concentrations of aerosols requires regular measurement of their optical properties such as aerosol optical depth (AOD).
Materials & Methods
Algeria is a large country with little knowledge of the spatial and temporal diversity of AOD, and the low spatial resolution of existing products makes it very difficult to predict aerosols (airborne particles) at the local scale, especially in arid southern regions. As a result, AOD recovery with data with higher spatial resolution is crucial for determining air pollution and air quality information. Several AERONET stations have been installed in Algeria. The Tamanrasset_INM station has been selected based on its location and the availability of historical AOD data for the period (2015-2016).
In this study, Landsat-8 / OLI image from tile 192/44 was used for satellite images. To this end, 23 TOA-corrected L1G-level Landsat-8 / OLI cloudless scenes were downloaded from January 2015 to December 2016 in the study area. DN values are converted to TOA reflections using the scaling factor coefficients in the OLI Landsat-8 metadata file. In this study, the minimum monthly reflectance technique was used to recover AOD in this area. As a result, LSR images were used in the recovery process in different months of 2015 and 2016. The process of selecting reference LSRs was initially based on the selection of clear, foggy / cloudless sky images. The selected images were then used to construct artificial images in which each pixel corresponds to the second lowest surface reflection of all selected monthly images to be the LSR pixel for the respective month. The AOD retrieval method developed in this study is based on a LUT, using the 6S radiative transfer model. The advantage of using the 6S model is its ability to estimate direct components and scattering using a limited number of inputs for each spectral band in the entire solar domain. The effect of the viewing angle is limited because Landsat data are usually obtained with a fixed viewing angle. Surface reflectance can be estimated from a pre-calculated LSR database. The accuracy of AOD recovery depends on the use of the appropriate aerosol model. A continental model was selected from the available aerosol models. Other atmospheric parameters such as ozone, carbon dioxide, carbon monoxide and water vapor are considered by default. The AOD values used to make LUT are set as follows: 0.0, 0.05, 0.1, 1.5, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.2 and 1.5. The zenith angles of the sun and the sensor range from 0 to 70 degrees with a step of 5 degrees and the range of azimuth angles from 0 to 180 degrees with a step of 12 degrees. Using these parameters, the radiative transfer equation was run in forward to obtain the TOA reflection. Different combinations of input and TOA output parameters are stored in LUT. AOD retrieval is based on a comparison between the TOAs estimated with the model and the observed items using the best fit approach. Using such an approach, the estimated AODs are simulated in accordance with those used in the production of TOAs, using a competency function that minimizes the distance.
Results & Discussion
In this study, the AODs recovered at 550 nm in a 5-by-5-pixel window around the AERONET site were averaged. The considered AERONET values are the average of all measurements taken within ± 30 minutes of image acquisition time. Observation regression results (AOD from Landsat 8 images and AERONET stations) showed that the correlation coefficient is about 84%. This study shows a good fit of the model on the research data and shows the high capability of the model. This study showed a strong recovery of AOD against AERONET data of more than 70% at . The differences can be attributed to a limited number of points or hypotheses related to the aerosol model used in this study. The assumption of using a pre-calculated LSR does not limit the accuracy of this method because we have shown that in arid regions where the change in land cover in different months of the year is small, a pre-calculated LSR image can be representation used the share of surface reflection in the radiative transfer model throughout the month.
Conclusion
In this study, an AOD derived from a high-resolution satellite at an urban scale was produced in the city of Tamanrasset, Algeria. The developed method assumes that the change in land cover is minimal and the temporal change in LSR is not significant. A pre-calculated LSR image is created to show the surface reflection in the retrieval process. Based on the 6S radiative transfer model, an LUT was constructed to simulate the TOA reflection of the built-in LSRs and a set of geometric and atmospheric parameters. The retrieved AODs were compared with the AERONET ground data. The results show that this approach can achieve reasonable accuracy in AOD recovery, which reaches about 70.9% at . In addition, this approach is suitable for estimating AOD in urban areas compared to existing AOD products with low spatial resolution. The results of this study show a 4% improvement compared to the results of Omari et al. (2019). The results of this study showed that ignoring the monthly changes in LSR values leads to good results in AOD recovery.
Saeed Farzaneh; Mohammad Ali Sharifi; Amir Abdolmaleki; Masood Dehvari
Abstract
Extended AbstractIntroductionSatellites in geodesy receive and transport important information. Among those, satellites with Low Earth Orbit (LEO), which are at altitudes less than 1000 km, have a significant role in the advancement of geophysical sciences such as earth’s potential field. ...
Read More
Extended AbstractIntroductionSatellites in geodesy receive and transport important information. Among those, satellites with Low Earth Orbit (LEO), which are at altitudes less than 1000 km, have a significant role in the advancement of geophysical sciences such as earth’s potential field. Many parameters have an impact on the precision and accuracy of their information. Atmospheric friction is one of the most principal forces on satellites, which may cause deviation and falling of satellite on a short period. From the beginning of aerospace missions, many efforts have been done to determine atmospheric friction by geodesists, e.g., empirical models of atmosphere neutral density. Because of the complex nature of atmosphere behavior and also data limitations, these models may have low accuracy. So, there is a need for methods to improve the accuracy of empirical models by means of combining observations of atmospheric density to predict its future state. Materials & MethodsAlong with the extension of computer science, new reliable algorithms have been introduced which are able to predict a time series; Artificial Intelligent (AI) and Neural Networks (NN) are the best of these methods. These simple algorithms are inspirations of the human brain and its ability to learn and have been used in many different scientific fields. In these techniques without any requirement for constructing complex modeling, the relation between input and output will be provided only using weight and bias vectors during the training procedure. Simple Neural Networks are memoryless meaning that the value of time-series in previous can’t be used for predicting the future value of time series and therefore some important dependency of signal values with time will be lost. A Recurrent Neural Network (RNN) has been implemented to overcome this issue. RNN’s can store some important information of the values of the time series in the previous steps in a chain-like structure and using this information for predicting the next value of time series that will improve the accuracy of prediction. In this study, the Long Short-Term Memory (LSTM) Neural Network which is a kind of Recurrent Neural Network’s has been implemented to predict the scale for correcting atmospheric density of numerical models. The data of Grace Accelerometer observation in the 6 first month of the year 2014 have been used for training the LSTM for univariate training. Also, the LSTM has been trained in multi-variants mode once with using the coefficient of atmospheric correction expansion up to degree 2 and once with using sun geomagnetic information along with information of k_p index. Results & DiscussionAfter training the LSTM network, by using the estimated parameters of the model, the zero degrees coefficient of harmonic expansion for a scale factor of correcting atmospheric density has been predicted in periods of 7, 14, 30, 60, and 90 days. The results of the univariate model show that the lower RMSE (Root Mean Square Error) is obtained about 0.054 in the period of prediction of about 14 days. Also, the results show that the multi-variants model with input data of sun geomagnetic information and k_p index has lower RMSE values in considered prediction periods compared to the other modes and the lowest RMSE is about 0.03 and belongs to the prediction of about 7 days. For evaluation of LSTM parameters in the obtained results, the predictions have been implemented with various Window sizes. The results show that by increasing windows size, the RMSE of the prediction will be reduced and the lowest RMSE was for prediction of 7 days with a window size of about 90 days. For the purpose of more evaluation, with the predicted atmospheric densities correction coefficient, the orbit of GRACE satellites has been propagated and the calculated position and velocity of satellites have been compared with the real orbit data. The results show that the lower RMSE will be provided with the prediction of 7 days with an RMSE for position and velocity of about 50 meters and 0.15 m/s respectively. ConclusionIn this study, due to the complex nature of the atmosphere, the LSTM Neural Network has been used for modeling and predict the zero-order scale for correcting atmospheric densities harmonic expansion. For training the network, the data of Grace Satellites Accelerometer in the 180 days of the year 2014 have been used. The LSTM has been in univariate and multi-variant models. In the multi-variants model, once with using the coefficient of atmospheric correction expansion up to degree two and once with using sun geomagnetic information along with information of k_p index the network have been trained. The period of prediction was considered of about 7, 14, 30, 60, and 90 days.The results show that the LSTM is capable to predict the correction coefficient in considered periods with a mean RMSE of about 0.05 for zero-order degree. Also, the results show that the lowest RMSE was for the 7 and 14 days of prediction and by increasing the window size of LSTM the RMSE will be decreased. The results of calculating the position of GRACE satellites position and velocity using predicted correction coefficients with real data show that the lowest RMSE was for prediction of 7 days for implemented method.
Saeed Farzaneh; Mohammad Ali Sharifi; Seyedeh Samira Talebi
Abstract
Extended Abstract
Introduction
In recent years, the development of the country in the space industry and the ability of building, launching and infusion of satellites in the lower orbit has put the limited number of countries with such technology. In order to complete the entire cycle of the space ...
Read More
Extended Abstract
Introduction
In recent years, the development of the country in the space industry and the ability of building, launching and infusion of satellites in the lower orbit has put the limited number of countries with such technology. In order to complete the entire cycle of the space industry, the satellite navigation and control, which has been neglected since the beginning of the movement of space science in the country, has to be considered specially. The attitude determination in one sentence is the application of a variety of techniques for estimating the attitude of spacecrafts. In dynamic astronomy, the attitude determination is the the process of controlling the orientation of an aerospace vehicle with respect to an inertial frame of reference or another entity such as the celestial sphere, certain fields, and nearby objects, etc.
A spacecraft attitude determination and control system typically uses a variety of sensors and actuators. Because attitude is described by three or more attitude variables, the di®erence between the desired and measured states is slightly more complicated than for a thermostat, or even for the position of the satellite in space. Furthermore, the mathematical analysis of attitude
determination is complicated by the fact that attitude determination is necessarily either underdetermined or overdetermined.
Materials and methods
Attitude determination typically requires finding three independent quantities, such as any minimal parameterization of the attitude matrix. The mathematics behind attitude determination can be broadly characterized into approaches that use stochastic analysis and approaches that do not. This paper considers a computationally efficient algorithm to optimally estimate the spacecraft attitude from vector observations taken at a single time, which is known as single-point or single-frame attitude determination method. There have been a number of attitude determination algorithms that compute optimal attitude of a spacecraft from various observation sources (known as the Wahba’s problem), and each of the methods has advantages and limitations in terms of accuracy and computational speed. The most popular are: the very important ˆq-Method, the most popular TRIAD and QUEST, SVD, FOAM, and ESOQ-1, the fastest ESOQ-2, and many others approaches introducing new insights or different characteristics, for instance, the EAA, Euler-2, Euler-ˆq, and OLAE.
Results and discussion
Since star detection algorithms can provide more than two stars, the star detector field of view often consists of two or more stars that are passed through the identification algorithms will be detected, those star vectors that have measurement errors can be compensated by using more than two stars. Methods such as the QUEST algorithm usually optimize an error function to the minimum optimal. In fact, the QUEST algorithm estimates the optimum specific eigenvalue and vector for the problem described in the Q_method method without the need for complex numerical calculations. The fact that the QUEST algorithm retains all the computational advantages of a fast definitive algorithm while maintaining the desired result efficiency underscores why it is typically used.
Conclusion
Simulation results showed that the traid and quest algorithms with shuster method attitude determination algorithm can be an efficient alternative over the eight tested algorithm in terms of computational efficiency for singularity-free attitude representation.
Saeed Farzaneh; Reza Shahhoseini; Iman Kordpour
Abstract
Introduction Drought is considered to be one of the most widespread natural disasters, ranking second in terms of damages. Due to the complex relationship between hydrological cycle parameters and atmospheric observations, predicting or modeling drought lacks the necessary precision. One of the most ...
Read More
Introduction Drought is considered to be one of the most widespread natural disasters, ranking second in terms of damages. Due to the complex relationship between hydrological cycle parameters and atmospheric observations, predicting or modeling drought lacks the necessary precision. One of the most significant problems in drought monitoring is lack of proper spatial coverage for the collected data (due to unavailibility of field data in some regions) and also lack of a suitable time scale (observations and thus drought estimation is not always possible). Since satellite observations do not face challenges like lack of spatial scale which is quite common in field observations, remote sensing satellites can provide a better estimate of droughts. However, satellite observations alone are not capable of accurately estimating the occurrence of droughts. Therefore, a combination of field and satellite observations has been used recentely to reach a better estimate of hydrological problems. Materials & Methods Temporal and spatial complexity of droughts have made a new global index combining ground-based and satellite-based observations quite necessary. Given the kind of data used in MDI index, we cannot expect it to be global. However, its performance is still acceptable in similar environments and climates, and thus it has been used in the United States (Texas). Datasets selected for the present study have different temporal and spatial scales and thus, a common scale must be found before calculating the index. Data received from GRACE satellite and MODIS sensor were downloaded monthly, but precipitation data were collected on a daily basis. Thus, aritmatic mean of precipitation data was calculated to reach a monthly avarage. Regarding the spatial scale, one-degree precipitation data were received from GRACE and MODIS while precipitation data extracted from synoptic stations had a point-based nature. Therefore, Inverse Distance Weighting (IDW) method was used to produce a one-degree network. Three types of observations were used in the present study including data received from synoptic stations of Iran meteorological organization, GRACE mission satellite-based gravity data and MODIS remote sensing satellite-based data. These were selected to identify droughts over a 14-year time series. Results & Discussion The present study has calculated MDI drought index on a one-degree spatial scale and monthly temporal scale for 168 months using Precipitation, NDVI, and TWS data. Severe droughts in northwestern and central areas of Iran from 2004 to 2014 have led to a shortage of water in reservoirs. In addition to drought, too much water harvesting in northwestern Iran has resulted in a decrease in groundwater level and thus, increased water harvesting from rivers and canals leading to the Urmia Lake and reduced water level in this lake. The results of MDI drought index calculated for Iran over the period of 2000 to 2014 show a high correlation with the results of standardized precipitation-evapotranspiration drought index. According to the type of data used to calculate MDI index, it is expected to have a strong correlation with PDSI index due to its sensitivity to precipitation, area temperature and soil moisture content. Since GRACE and MODIS satellite-based data, and data received from synoptic stations were used, a strong correlation with MDI is also expected. It should be noted that PDSI index is higher than MDI index in Iran, although both show the drought trends accurately. For example according to PDSI index, the worst drought of the last two decades in Iran has occurred in 2008, and MDI index shows the same year. Conclusion The present study has introduced a new drought index using a combination of precipitation data, GRACE_TWS and NDVI. These data were selected because of their high sensitivity to drought. GRACE_TWS observations monitor hydrological drought and include surface and subsurface water sources. NDVI observations are mostly used to identify photosynthetic activities of vegetation cover and are therefore very useful for detecting agricultural drought. Precipitation value shows the amount of surface water in the study area. Precipitation can have relatively rapid effects and is therefore useful for monitoring meteorological drought. MDI index has identified several droughts in each region of the country in the period of 2003 to 2016. These identified droughts have generally covered the country over time. However, each drought has had a different impact on ecosystem. In Iran, the most severe droughts have occurred during 2008 to 2009 and 2011 to 2012. Since MDI correlates well with PDSI, both show a drought in these years. In order to develop the proposed algorithm, the effect of different zoning of the study area on MDI index can be studied.
Zahra Banimostafavi; Saeed Farzaneh; Mohammad Ali Sharifi
Abstract
Extended Abstract:
Introduction
Nowadays, engineering structures face many threats. Natural and human activities can result in deformation and displacement of dams, bridges, and towers. As a result, any crack in the body of these structures is important and may have dangerous consequences. To prevent ...
Read More
Extended Abstract:
Introduction
Nowadays, engineering structures face many threats. Natural and human activities can result in deformation and displacement of dams, bridges, and towers. As a result, any crack in the body of these structures is important and may have dangerous consequences. To prevent catastrophes, the behavior of these structures should be monitored permanently during the construction phase and after opening.Nowadays,thebehavior of engineering structures such as dams, power plants, and towers is considered to be especially important. Three different methods are usually used to measure such behavior: classical, satellite and precise instruments.
Materials and methods
Modern equipment is considered to be a crucial factor in controllingpossible changes and preventing human errors. Therefore, different sensors are installed in the structure to measure tensile and shear flexibility during the construction phase. Moreover, data received from these sensors is analyzed permanently during the service life to ensure sustainability of the structure. These tools make internal analysis of these structures possible. Analyzing the behavior of engineering structures is considered to be one of the most important tasks in the field of geodesy. Inaccurate analysis of displacements can have deadly effects. Various methods are used to measure such displacements, which are divided into two categories: robust and non-robust methods based upon the results of the epoch adjustment. To find deformations, a geodetic network should be defined in the first step. If two epochs are not measured in the same datum, the results will not be reliable. Displacement can be measured in two ways: Absolute and Relative. In the absolute method, some points are considered to be stable, while in the relative network, all points are considered to be unstable, and the problem is solved based upon this hypothesis. The method of relative network is used in the present study. Regarding network geometry, displacement analysis is performed using two methods:single and combinatorial. Moreover, displacement analysis is divided into two categories of robust and non-robust methods. Iterative Weighted Similarity Transformation (IWST)and Minimum L1 norm are among robust methods which calculate the matrix of displacement by minimizing the first and second norm. Global Congruency Test (GCT) is a non-robust statistical method used to determine unstable points in geodetic networks. Robust and GCT are among classical methods used to discover unstable points in geodetic networks, while Simultaneous Adjustment of Two Epoch (SATE(is a new method used to achieve this purpose. Combinatorial methods are also considered to be a suitable alternative method used for detecting unstable points in a geodetic network. In our previous study, “evaluation of single-point methods used fordetecting displacement in classical geodetic networks”, single-point methods of detecting unstable points were investigated and the SATE method was selected as the optimal method. Unlike single-point methods, these methods examine all points of the geodetic network simultaneously to discover unstable points.
Results and discussion
The strong dependence of these methods on the network geometry makes discovery of all unstable points impossible. Combinatorial methods are considered to be a suitable alternative method used to detect all unstable points in the geodetic network. These methods does not have a strong dependence on scale and the network geometry. Multiple Sub Sample and M-split methods are classified in this category. These methods can detect unstable points efficiently. The present study takes advantage of simulated datato evaluate combinatorial methods such as Multiple Sub Sample (MSS) Angles, MSS-distance difference, and M-split and compare them with the SATE method with the aim of choosing the optimal method. Then, unstable points in the real network of Jamishan dam in Kermanshah Province will be discovered using the identified optimal method.
Conclusion
The present study identifies the best method between single and combinatorial methods. The best method can detect most unstable points and has the lowest dependence on geometry, scale and other factors influencing the results.According to the results, Multiple Sub Sample with distance difference is selected as the best method.